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1. INTRODUCTION 

Ideal power system network is expected to operate with sinusoidal voltages and currents of a single 
frequency and constant amplitude. This ideal operation is not possible in realistic power systems because of 
increasing use of power electronics devices in modern power system network. Nonlinear components and 
loads, such as power converters, Flexible AC Transmission Systems (FACTS) devices, nonlinear saturation 
and hysteresis in components with magnetic cores represent the main contribution harmonic distortion. 
Flexible alternating current transmission systems devices are widely used in transmission network to improve 
electrical and economic performance. But they are nonlinear devices, they generate undesirable distortions in 
voltage or current waveforms. Distorted waveforms contain unwanted harmonics. This harmonic distortion 
creates many problems such as misoperation of protective devices, reduces the life of electrical equipment, 
and increases losses in the power systems [1]. Considerable progress has been made in the development of 
correct apparatus to monitor the harmonic behaviour of power systems [2, 3]. However, in system planning 
and system analysis problem must be considered separately because measurements may not be economic or 
the system may not be available. In such situations, computer simulations based on the mathematical 
modelling provide the workable alternative to physical measurements [1]. Many works have been proposed 
mathematical models of nonlinear devices to determine waveform distortion in power systems. To study 
the steady-state voltage and current waveform distortion, time domain simulation has been developed many 
years back. This technique usually takes long simulation run times to calculate the harmonic content of 
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distorted waveforms even for small systems [4]. This is because time domain technique first finds steady 
state condition after that harmonic content of distorted waveform is found using the Fast Fourier Transform 
(FFT). The complications in time domain provided motivation for developing the frequency domain 
techniques. Frequency domain techniques take less simulation time as compare to the time domain 
techniques. The most fruitful frequency domain technique is Harmonic Domain (HD) technique which is 
used for both analysis and modelling of harmonic sources. Complex Fourier Series (CFS) and Norton 
equivalent circuit formulations are used in Harmonic domain technique. HD solution process have presented 
in reference [1]. 

HD has been applied to obtain steady state response of the power system components such as 
synchronous machines, the power transformer, arc furnaces, Fluorescent lamps, Thyristor Commutated 
Reactors (TCRs) .The Harmonic Domain Technique has been applied to model the Thyristor Controlled 
Series Capacitor (TCSC) for both stability analysis and steady state analysis. In this paper, Harmonic Domain 
approach is used for the harmonic analysis of three phase delta connected Fixed Capacitor-Thyristor 
Controlled Reactor operating under balanced and unbalanced conditions. The paper is organized as follows. 
Basic operation of Single phase FC-TCR is presented in the second section. The third section describes 
modelling of three phase FC-TCR in HD approach. Test system description is given in the fourth section. 
Simulation results are presented in the fifth section. 


2. BASIC OPERATION OF SINGLE PHASE FC-TCR 

This section gives the main components of FC-TCR and basic operation of single phase FC-TCR. Figure 1 
shows the single phase FC-TCR used to supply variable reactive power support at a specific location in 
Power systems. Basically, single FC-TCR consists of fixed capacitor, reactor and two thyristors. Connections 
of basic elements are shown in Figure 1. It is possible to operate the FC-TCR in both capacitive mode and 
inductive mode with help of firing signals given to thyristors. Firing angle a plays important role to control 
reactive power. Firing angle æ is varied from 90° to 180°. TCR is full conducting state at a=90°. TCR is said 
to be full conducting mode when firing angle (a) is 90° Figure 2 shows the TCR current (Iter) and 
TCRvoltage (Vicr) waveforms when the TCR is in full conducting mode. 
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Figure 1. Single phase FC-TCR 
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Figure 2. TCR Voltage and TCR current waveforms when a@=90° 
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It is clear from Figure 2 that there is no distortion in current waveform. No current harmonics are 
generated when TCR in full conduction mode. TCR is said tobe non-conducting mode when firing angle 
(a) is 180°. Figure 3 shows the TCR current and TCR voltage waveforms when TCR is non-conducting 
mode. It is clear from Figure 3 that there is no current flowing in TCR branch. Because both thyristers are 
non conducting. TCR is said to be partially conducting mode when firing angle is varied between 90° to 180° 
Figure 4 shows the TCR current and TCR voltage waveforms when TCR is partially conducting mode. 
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Figure 3. TCR Voltage and TCR current waveforms when a=180° 
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Figure 4. TCR Voltage and TCR current waveforms when a=120° 


It is clear from Figure 4 that there is distortion in current waveform. When TCR is partially 
conducting, it will generate the harmonics at connection point. One way to obtain harmonic content 
generated by FC-TCR is measurement by power quality analyser. We can predict harmonic content of 
FC-TCR using computer simulation.Simulation techniques does not requires the costly power quality 
analysers. 


3. MODELLING OF THREE PHASE FC-TCR 

Peter wood proposed the switching function concept to determine the harmonic content of 
converters. Later Bohmann and Lasseter the same switching function method applied to basic TCR circuit to 
compute the harmonic content of TCR. The main drawback of this method is that switching instants are 
fixed. Acha introduced the Norton equivalent circuit model of TCR. Harmonic content of switching function 
depends on conduction of thyristors. The switching function of TCR is a function of TCR voltage and firing 
anglea. The TCR voltage zero crossing is appropriated as the reference for giving the firing signal. 
Mathematical modelling of TCR should include the thyristors turn-on and turn-off instants. Harmonic 
domain model of single phase FC-TCR is given by the following Equation (1) [5, 6]. 
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Yrcr is admittance matrix of TCR which represents the steady state operation of TCR as shown in 
Equation (2) and (3). 


Yrer = +D- (jhw)S (2) 
Yc = CD(jhw) (3) 


S is harmonic switching vector. D is differention operator. Practical power systems, three phase FC-TCR 
consists of three single-phase TCR are connected in delta and three capacitors are connected in star. 
The equations for three phaseTCR is given by Equation (4). 











Ir F Yrer + Yfir —Yrer Yrer] [VR 
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4. SIMULATION RESULTS 
Three phase FC-TCR model developed in section 3 is tested on two bus system shown in Figure 5. 
The system contains the FC-TCR at bus 2. Figure 5 is reduced to simple circuit as shown in Figure 6. 
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Figure 5. Test system under consideration Figure 6. Equivalent circuit of two bus system 























V, is source voltage. V2 is voltage at bus2.Three phase FC-TCR is connected at bus 2. V2 is unknown which 
can be deternine form network equations. Network Equations can be written for Figure 6 as follows. 


Network Equation (5) 
Pee Ye22 | k] (5) 
0 Ye21 Y:22 + Yrc-rcr] LVzryB 
Virypg is source voltage. V2pyp Can be determine from Equation (6) iteratively 
Vzryg = —(¥t22 + Yrc-rcr) 'Yt,22 Virve (6) 


The following parameters are used to simulation. The Transformer parameters: Star/delta, 
220kv/33kv, 10% of impedance, 100 MVA. FC-TCR rating: (0-80) MVAR. Three cases are considered here 
and simulation results have been presented. 
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4.1. Balanced operation of FC-TCR 
Three phase voltage magnitude are equal and phase displacement is 120°. Results are presented for 
two firing angles under this case. 


4.1.1. Firing angle=95° 

Figure 7 shows tha harmonic content of voltages at bus 2.Magnitude of Fifth harmonic is more as 
compared to the other harnonics. This because of TCR current waveform, 5" harmonic is dominant which is 
clear from the Figure 8. Fifth harmonics is dominating in capacitor current waveform as copmpare to 
the other harmonics. 


4.1.2. Firing angle=100° 

Fifth harmonic is dominating here also magnitude is more as compare to previous case. 
This because of TCR current waveform becomes more nonsinusoidal as firing angle increases. It is also 
observe from Figure 7 to Figure 12 that 3 harmonic is not appeared because 3" harmonic produced in TCR 
is circulated in delta connected TCR under balanced condition. 
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Figure 7. Three phase voltages harmonic spectrum at bus 2 
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Figure 8. Three phase TCR currents harmonic spectrum 
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Figure 9. Three phase Capacitor currents harmonic spectrum 
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Figure 10. Three phase voltages harmonic spectrum at bus 2 
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Figure 11. Three phase TCR currents harmonic spectrum 
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Figure 12. Three phase Capacitor currents harmonic spectrum 


4.2. Unbalance in voltage magnitudes 

Unbalance in voltage magnitudes assume voltage magnitude values in three phase voltages as 
Vr=1 P.U, Vy=0.95pu, Vp=1.2pu values and Firing angle=100° are considered. Harmonic content of voltages 
at bus2, TCR currents, capacitor currents are shown in Figure 13, Figure 14 and Figure 15 respectively. 
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Figure 13. Three phase voltages harmonic spectrum at bus 2 
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Figure 14. Three phase TCR currents harmonic spectrum 
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Figure 15. Three phase Capacitor currents harmonic spectrum 


It is clear that 3" harmonic is appeared when there is unbalance due to voltage magnitudes and 
magnitude of 3"! harmonic is less as compared to the 5" harmonic magnitude. 


4.3. Unbalance in firing angles 

Unbalance in firing angles votage magnitude of all three phase voltages are equal and different 
firing angles are considered in this case. Firing angle in R phase =95°, Firing angle in Y phase =100°, Firing 
angle in B phase =120°, and equal voltage magnitudes (1pu) are used in simulation. Results are presented in 
Figure 16, Figure 17 and Figure 18. 








E E E E 
HEE F phase 
[Ay phase 
2.5 HE £ Phase H 























magnitude 
a 














0.5 


Hu 


alle Pant 
2 3 4 5 6 7 8 9 10 11 12 13 14 15 
harmonic order 












































Figure 16. Three phase voltages harmonic spectrum at bus 2 
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Figure 17. Three phase TCR currents harmonic spectrum 
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Figure 18. Three phase Capacitor currents harmonic spectrum 


It is clear that 3"? harmonic is also appeared here also and magnitude of 3 harmonic is more as 
compared to the 5“ harmonic magnitude. 


5. CONCLUSION 

Harmonic analysis of three-phase fixed capacitor-Thyristor Controlled Reactor working under 
balanced and unbalanced conditions has been investigated using HD approach. The results shows that 
significant amount of 5" harmonic appears in voltage waveform at connection point of FC-TCR, TCR per 
phase current and capacitor line current when system is operating under balanced conditions. 3 harmonic is 
not appeared when system is operating under balanced condition. Third harmonic is appeared when there is 
unbalance due to voltage magnitudes and magnitude of 3"! harmonic is less as compared to the 5" harmonic 
magnitude. Third harmonic is also appeared when there is unbalance due to firing angles and magnitude of 
3" harmonic is more as compared to the 5“ harmonic magnitude. This work is useful to design the filters to 
mitigate the harmonics under different operating conditions. Steady state harmonic analysis is carried out in 
this paper. Future scope of this work is to determine harmonics of given system under transient conditions. 
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